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In this letter, it is pointed out that the two matrix model defined by the action S = 
\ (tr + tr S^) - ^ tr - ^ tr S"^ - I tr( AS) ^ can be solved in the large limit using 
a generalization of the solution of Kazakov and Zinn- Justin (who considered the symmetric 
case a A = cts)- This model could have useful applications to 3D Lorentzian gravity. 
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1. Introduction 

There are a variety of motivations to study the large N limit of matrix integrals, from 
statistical mechanics on random lattices and graph enumeration [0 to quantum chaos and 
mesoscopic systems 0. In [Q], the so-called ABAB model was introduced, as an example 
of model exhibiting the critical behavior of c = 1 conformal field theory coupled to gravity. 
It has the partition function 



Z(a,/3) = dAdB expN 



- - (tr ^2 + tr S^) + ^ (tr + tr S^) + - tr{ABf 
2 2 



(1.1) 



The special form of the interaction between the two matrices, tr ABAB, which is different 
from the usual two-matrix model Q , makes it impossible to compute this partition function 
using the standard method - direct diagonalization of A and B. Instead a new method was 
proposed, combining several techniques which had appeared in previous articles, including 
character expansion [|5|, saddle point over Young diagrams and functional inversion 
relations 0. It led to an exact expression for the "planar" (large A^) free energy in terms 
of elliptic functions. Subsequently, this model was used in several contexts: 2D Eulerian 
gravity , knot theory and even 3D Lorentzian gravity [|l^ . This prompted the author 
to consider the following natural generalization: 



Z(a^,aB,/3)= dAdBexpN 



--{trA^ + trB') + ^ tr + ^ tr + ^ tr {ABf 
2 ^4 4 2 

(1.2) 

in which the coefficients of the quartic terms are completely arbitrary. The Feynman rules 
of this model are as follows: we form fat graphs with edges in two colors (red and green) 
and quartic vertices, in such a way that colors can only cross (each other, or themselves) 
at each vertex. Red (resp. green, mixed) vertices are assigned a weight of a a (resp. a^, 
/3). Alternatively, let us note that if one sets X = then one obtains 



Z{b,c,d,e) = J dXdX^ expN 



- tr{XX^) +btr{X^X^^) + - tr{XX^)^ 
+ ^ tv{X^ + X^^) + etv{XX^^ + X^X^) 



(1.3) 



with b = (aA + o:B + '2l3)/4, c = d = {aA + ctB — 2/3)/4 and e = (a^ — q;b)/4. It corresponds 
to a perturbation of the symmetric case a a = olb (which describes an 8- vertex model) with 
the extra operators X'^X'^ -f h.c, turning it into a 16- vertex model on random dynamical 
lattices. We shall come back in the concluding section to possible applications of this 
model. For now, let us show how, despite the more complicated analytic structure than in 
the symmetric case, one can solve it. 
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2. Character expansion and saddle point equations 



As in 0, we immediately transform the expression ( p. .21) by means of a character 
expansion of the term exp[A^^ tr{AB)^] as a function of AB. We obtain: 



7r„ r, !J) - Vr.V.q.'^l*'-/^ A(ft"«V2) A((h-^-^ - l)/2) 



dA dB exp 



- (tr A2 + tr + ^ tr + ^ tr B' 
2^ ^4 4 



f2.1) 



XhiAB) 



where the sum is over a set h 
0, and = ^hi— 



{hi, . . . , h]^} of integers that satisfy /ii > /i2 > • • • > /^at > 
is the number of boxes of the Young diagram. A(-) is the Van 



der Monde determinant, Xh is the GL{N) character associated to the set of shifted highest 
weights /i, and the f^:^^^/"'^'^ qj^-q ^j^g even/odd hi, which must be in equal numbers. It is 
now possible, using character orthogonality relations, to integrate over the relative angle 
between A and B; this leads to a separation into one-matrix integrals: 



Z(a^,aB,/3) = Y,iN(3/2)*''/''chRhicxA)Rh{aB) 



(2.2) 



where Ch is a coefficient: Ch 
integral 

Rh(a) = 



dMxh{M)expN 



7 ^ 



and Rhid) is the one-matrix 



--trM^ + -trM^ 
2 4 



(2.3) 



which was studied in detail in [Q. 

We can now consider the large N limit: the summation over the hj is dominated 
by a saddle point which is characterized by a continuous density p{h) of the rescaled 
variables hj/N. There is generically a saturated region [0, where p{h) = 1 and an 
unsaturated region where < p{h) < 1. In order to write down the saddle point 

equation we need to differentiate with respect to the hi] we denote by d /dh the operation 
(l/A)a/a/ii inthelarge Alimitwith/ii/A = h. Then we have 2(a/a/i) logi?/^(ayi) =^(/i), 
2{d/dh) \ogRh{cxB) = ]^{h), (d/dh) \og{ch/N*^/'^) = -{${h) + log |)/2, where the slash 
notation means Ip{h) = ^{D{h + iO) +D{h — iO)) for any analytic function D{h) with a cut 
on the real axis, L{h) and M{h) are unknown functions whose analytic structure will be 
discussed below, and H{h) = limAr^oo J2i h-h 
saddle point equation: 



/ '^'h-h' • This results in the following 



l^{h) +I^{h)-${h) 



log^ /lG[/li,/l2] 



(2.4) 
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3. Analytic structure and solution 

From the study of Rh{ot) performed in and that we shall not repeat here, it follows 
that the function L{h) defined by 2{d / dh) log Rh{(y- a) = Jl'ih) on [/ii,/i2] and by the fact 
that it has the same cut as H{h) on [hi, /i2], has the following analytic structure: besides 
a logarithmic cut starting at /i = and the cut of H{h) on [hi, /i2] (i.e. discontinuity of 
2nip{h)), it has a semi-infinite square root cut starting at h = hA, which we choose to be 
[hA, +oo]. Similary, M{h) defined by replacing with as has a cut starting ai h = hs- 
Furthermore, they satisfy the following equations: 

21^{h) = H{h) + log— he[hA:+oo] (3.1a) 
aA 

21^ (h) = H{h) + log— hE[hB,+oo] (3.16) 

as 

Now define the function F{h) = L{h) + M{h) — H{h) — log h. This is a rather compli- 
cated function since it has 3 square root cuts leading to other sheets. The equations ( pT^ ) 
and ( pT|) define the monodromy of F{h) around the branch points h±, hA and hs- If one 
differentiates once F{h), the inhomogeneities in Eqs. ( p.4|) and (|3.1j ) disappear (as well as 
the logarithmic cut) and we find that fortunately, F'{h) has only 4 sheets: the original 
sheet F'{h), its opposite —F'{h) by going through [hi,h2], and ±G"(/i) by going around 
hA or hs, with G{h) = L{h) — M{h). Therefore it is a meromorphic function of u, where 

h-^-^±^+u^+^^^^^^=Q (3.2) 

This shows that the problem is solvable. In practice, to stay closer to the solution of 0, 
we introduce z = —v? = ^{h — ^^^^^ ± {h — hA){h — hs)) , with a cut on [hA, hs], the 
-f sign being for the sheet of ±F, the — sign being for ±G; also introduce the images z±, 
zq, za, zb on the F sheet of h±, 0, hA, hs- Then the function F{z) (which combines F{h) 
and G{h)) has square root cuts on z+] and [0, +oo] and satisfies 

^{z) = -log(3 ze [z- , z+] (3.3a) 
^{z) = -^logiaAas) z G [0, +oo] (3.36) 

Taking into account the extra logarithmic cut on [zq, we find that F{z) is of the form: 
F(z) = +2*,,(.) - (3.4) 



where $5(2;) is the incomplete eUiptic integral of the third kind, which can be expressed 
explicitly as follows: define /c^ = the complete elliptic integrals K, K' , E, and 

the parameterization sn^(x, k) = z-j z. Then ^i{z) = x{2Zi{x) — ^) +log ^^jf^fj, where 
Oa (a = 1, 2, 3, 4) is a rescaled theta function and Za its logarithmic derivative. 

The parameters h±, Ha, hs are fixed by computing the asymptotic behavior of F{z) 
as z 0, z ^ 00 and the normalization of the density J dhp{h). This calculation is 
similar (though more complicated) than what is performed in and we skip the details. 

First, based on the known asymptotic behavior L{h) = ^ log ^ — 2y/-aAh 2k 
0(/i"^/2) (and similarly for M{h)) we obtain 3 conditions as z ^ 00, 0: 

^0 = y + K (3.5a) 

n^ = ^-^{a-//' + a-'/') (3.56) 
-1/2 _ -1/2 

O4 = 2 'y (3.5c) 

'I' A — il^B 



where 

log(/32/aAaB) , .7/ N 
= + 4Z„(xo) (3.6) 

and xo is defined by sn^(xo, k) = Z-/ zq. 

The normalization of the density, / dhp{h) = 1, after integration by parts, is essen- 
tially a complete elliptic integral of the second kind: 

zoK'n^- z-E'ni + ^K'fli- ^—fl4 = 7i (3.7) 

Zq Z- k^ 

Equations (p.5|)-(|37^) provide an explicit parameterization of the quantities of the 
model in terms of k and ratios of coupling constants, such as (3'^ /aAdB and aA/cts- For 
example, we find 



n\ 40f V V"^/ 40l V V"^ 



(3.8) 

where we recall that xq and Qa are known via Eqs. (|3.5a|) and (p^). 

From the function F{z) one can derive many correlation functions, including 
(|^(Ai?)^) and (^A'^), but the formulae are quite cumbersome and are not reproduced 
here. 
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4. Phase diagram 

Ordinary criticality, wtiicti in the context of random dynamical graphs means the 
emergence of very large graphs i.e. the continuum limit, appears when one of the square 
root singularities at /i = /i^, h = hs-, h = h+ degenerate into power 3/2 singularities. 
Thus, there are three phases denoted by lA, IB and II, which all correspond to pure 
gravity (c = 0). Criticality can easily be tested numerically and the result is Fig. 1. On 
this two-dimensional critical surface, there are critical lines corresponding to the transition 
between two of these phases and where the corresponding branch points coalesce; and these 
three lines meet at a critical point where all three branch points merge. 
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Fig. 1: Critical surface in the {aA, ob, (3) space. 



The easiest transition to analyze is lA/IB: Ha = hs implies a a = c^b and we go 
back to the symmetric model. According to the analysis of 0], at this transition line the 
model is still pure gravity. More interesting is the lA/II transition. It is given by the 
condition Ha = h^: note that contrary to the symmetric case, Ha = h+ does not imply 
that /c — i> (trigonometric limit of the elliptic functions). However this condition still 
produces a degenerate analytic structure where two square root branch points turn into a 
cubic branch point: this is typical of a c = 1/2 (i.e. critical Ising) theory. This is confirmed 



by the analysis of the vicinity of the critical point that comes now (also, based on it 
is reasonable to assume that the particular point as = of this line is the critical point 
of a modified 0(1) model, which is known to have c = 1/2). 

Indeed, the critical point a a = ct_B = /3 = 4^ was already studied in ||^. It describes a 
c = 1 model - a free boson coupled to gravity. However the vicinity, in the 3-dimensional 
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T-A = -^i^rrih^'^^ - + 2?)') + ■ ■ • (4-1) 



parameter space, displays a rather complex behavior since 3 critical lines merge there. 
Let us briefly analyze it. We deflne the following deviation from criticality: A = 1 — 

-i 1/2 1/2 

(a^ + ) . The critical point corresponds to the trigonometric limit where the 
elliptic nome q goes to zero. By small q expansion of Eqs. (|3.5|) -(|377|), we can express A 
as a function of q and of s = s^^"^'' ^ ~ (Trctyi)"^^^ — (Trcts)"^/^; differentiating, 
we obtain 

d ^ log q s — 2q 
dq ~ 32q {s + 2qY 

This shows the following scaling behaviors: s oc q and t oc Q"^/^, as well as defines the 
phases in the vicinity of the critical point: s = 2q for phase II, t = ±8q'~^/^(s + 2q)'^ for 
phases lA/IB. The line separating lA and II is therefore given by s = 2q, t = 128q'^^'^. 
Computing A on the critical surface allows to conclude that transition lA/IB is first order, 
whereas transitions lA/II and IB/II are third order. 

Finally, a typical correlation function F follows the same pattern as A itself. In 
particular on a generic point of the critical surface, ^F = and therefore F = A'^/^ + 
regular, which confirms that the central charge c is zero. However at the line separating lA 
and II one easily finds = = and therefore F = A'^/"^ + regular which according 

to the KPZ relation [|l^ means that c = 1/2. 



5. Concluding comments 

The asymmetric ABAB model which we have just solved is enlightening in several 
respects. First it is an interesting statistical model with a rich phase diagram which might 
deserve further investigation. But it is also, at a technical level, the first solution of a matrix 
model in which we have a non-trivial three cut structure (the multi-cut phases of standard 
matrix models usually have a very simple - e.g. hyperelliptic - structure) . It can be hoped 
that for many other matrix models, the seemingly complicated monodromy equations close 
and allow to provide an exact solution. Finally, at the physical level we expect that the 
asymmetric ABAB model may be of relevance to three-dimensional Lorentzian gravity. 
Indeed, it is precisely the model that was introduced in Jl^ (cf Eq. (23)) to describe time 
slices of three-dimensional geometries. We hope that the exact solution can shed some 
light on some of the remaining unsolved issues. 
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